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Based only on simple principles of renormalization in coordinate space, we derive closed renor- 
malized amplitudes and renormalization group constants at 1- and 2-loop orders for scalar field 
theories in general backgrounds. This is achieved through a generic renormalization procedure we 
develop exploiting the central idea behind differential renormalization, which needs as only inputs 
the propagator and the appropriate laplacian for the backgrounds in question. We work out this 
generic coordinate space renormalization in some detail, and subsequently back it up with specific 
calculations for scalar theories both on curved backgrounds, manifestly preserving diffeomorphism 
, invariance, and at finite temperature. 
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I. GENERIC RENORMALIZATION OF SHORT-DISTANCE SINGULARITIES 



>■ - 

Quantum field theories are often considered under classical external backgrounds. Two standard examples are 
thermal baths jjj and curved space-times |2j and, in general, the procedure of renormalization becomes much more 
involved due to the appearance of new dimcnsionful scales in the problem. 

It is nevertheless known that some of the first coefficients of the renormalization group constants remain unaltered. 
For instance, finite temperature can be encoded by compactifying the Euclidean time, and this obviously affects the 
long-distance properties of the theory but not the leading short-distance ones. It is natural to expect no modification 
of the first coefficient of the /3-function. Similarly, gravitational backgrounds do not affect the leading renormalization 
of flat field theories^] due to the Equivalence Principle, which allows for a locally flat space-time and thus guarantees 
that short-distance singularities are those of Minkowski space-time. On the other hand, the complicated structure 
Oh, of thermal and gravitational backgrounds makes it often impossible in practice to get closed expressions for even 
the simplest Green functions. Infrared modifications of a theory do not interfere with the renormalization of leading 
singularities but make almost impractical any computation at and beyond one loop. 

It is our aim to present a renormalization procedure adapted to quantum field theory defined at finite temperature 
or in curved space-times which produces closed explicit expressions for Green functions at low orders of perturbation 
theory. The method is based on the ideas behind differential renormalization (DR hereafter) ||, and we shall now 
expose it in a general setting. By way of example, we shall also check that the first two coefficients of the /3-function 
and the first one of the anomalous dimension of the scalar field in A</> 4 theory are independent of long-distance aspects 
of the theory. 

Let us consider a massless Euclidean scalar field theory with a propagator given by 

< 0(x)0(O) > ee G(x,0;a) , (1.1) 

where a stands for dimcnsionful parameters such as temperature or curvature. The propagator obeys an equation of 
the type 

VG{x,Q-a) = ^=±, (1.2) 
V9W 
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where T> is the proper scalar laplacian of the theory, for instance on curved space-time (with background metric 
9fj,v( x )) or with nontrivial boundary conditions when necessary. For a preview of a particular form of this operator, 
one may consider for example 

V = -gTix^^du - a 2 , (1.3) 

where V M stands for a covariant derivative and a 2 is related to the curvature of the background metric (see Sec. 2). 
For the sake of simplicity, we shall first consider massless A0 4 theory. The massive case is discussed at the end of this 
section. We shall assume that our field theory is locally flat, that is, 

G(x,0;a) ^ -LI (1.4) 

This is the case of field theories both in curved spaces or in thermal baths. The long-distance behavior of the 
propagator will surely be complicated. This is at the heart of the problem of defining a Fourier transform and 
computing it. 




Figure 1. 

As we turn to the simplest quantum correction of the amputated four-point vertex function, Fig. 1, we find the 
contribution 

y {G{ Xl Q-a)) 2 , (1.5) 
which displays a leading singularity at short distances of the type 

(G(x,0;a)f ^ (L6) 



Due to this (logarithmic) singularity at x = 0, Eq. (1.5) does not accept a Fourier transform, it is not a good distri- 
bution upon integration against plane waves and needs to be renormalized. Typically, the short-distance expansion in 
Eq. ( |1.6D will contain subleading terms which diverge as the propagator itself and need not undergo renormalization. 
In the absence of nontrivial backgrou nds , the method of differential renormalization gives a recipe to obtain right 



away the renormalized version of Eq. (l.E). The idea is to rewrite this equation by extracting a laplacian operator, 



1 \nx 2 M 2 , s 

lU—t—, (1.7) 



(a; 2 ) 2 4 x 

where □ = d^d^. This expression is an identity at x ^ and, furthermore, produces an extension of the too singular 
function of the l.h.s. into a proper distribution on the r.h.s., provided the operator □ is understood as acting onto 
the left. We should now generalize this simple identity to more complicated cases, when the infrared behavior of the 
propagator becomes very involved. We achieve this goal by noting that, in general, the logarithmic singularity in Eq. 



(1.5) is solved by writing it as 



G 2 (x) = --Ld (G(x) \nG(x)/M 2 ) + D(x) . (1.8) 

l07T z 

D(x) = D(x,0;a) is then a bona fide distribution which will become explicit in each particular case, and G(x) is 
short-hand for G(x, 0; a). Again, the function on the l.h.s. is turned into a distribution while nothing is changed away 
from the singularity. It is easy to see, by using the leading beha vior of the propagator, that indeed the above formula 
encodes the renormalization of the singularity as shown in Eq. fll.7| ). The form of D(x) will depend on the particular 
problem we are addressing, and may look involved since it is encoding long-distance properties of the theory, but it 
is guaranteed to be at most as singular as the propagator itself. 
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There are basically three reasons why Eq. ([l]^) produces the correct renormalization of the 1-loop diagram in 
thermal or gravitational deformations of massless \<j> theory: 

1) Removal of singularities. As we just explained, our generalized re nor malization of the bubble diagram is guar- 
anteed by the short-distance limit of all the pieces of the r.h.s. of Eq. (1.8). 

2) Renormalization group behavior. The renormalization scale only appears within the logarithm. The renormal- 
ization group equation at one loop^j 



M ^ +/?(A) J^ (A) ) r(4) U - 



(1.9) 



where 7(A) stands for the anomalous dimension of the 
four-point function, will thus contain the piece 



field, /3(A) is the (3- function and I^ 4 ) is the amputated 



d 
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M dM I "1*3 In G(z)/M 2 ) + D(z) ) = — 2 VG{x) = 



16tt 



8tt 2 



1 



8tt- 



(1.10) 



The full diagram carries an extra A 2 /2 factor plus the addition of the s, t and u channels. Therefore, the M dependence 
can be reabsorbed by a change of the tree-level coupling dictated by the beta function 



/3(A) = M 



OX 



3A 2 



dM 16tt 2 ' 



(1.11) 



and the anomalous dimension receives no contribution from this diagram. 

3) Unitarity. T he im aginary part of the amplitude is related to the cross-section through the standard unitarity 
relations. In Eq. ( jl.8| ), the imaginary part is carried by the logarithm, which in turn gives rise to a delta function. 
As expected, the imaginary part of the 1-loop contribution is then proportional to the tree level structure. 

It is also noteworthy that in the presence of a gravitational background, our procedure preserves diffeomorphism 
invariance since only covariant derivatives are manipulated. There is just the desired scale invariance breakdown. On 
the other hand, the above recipe has to be slightly modified in the case of massive theories as we discuss at the end 
of this section. 




Figure 2. 




Figure 3. 



At two loops, the amputated 4-point function is given by the two diagrams in Figs. 2,3. The first diagram 
corresponds to a convolution and, thus, introduces no new kind of singularity. It just produces a promotion of leading 
logs to second order. To see this, note that the bare amplitude 



2 New /3-mnctions associated to mass and background couplings will eventually appear, but at higher loop order. 



3 



d A w y/g{w) (G(x, w; a)) 2 (G(w, 0; a)) 



(1.12) 



simply gets renormalized using Eq. (1.8). The fact that this diagram brings no contribution to the 2-loop (3- function 
stems from the fact that its renormalization scale dependence is non-local, which is readily checked 



M 



a 

dM 



d^w Jg [ V( G( x, w; a) In G(x,w; a) /M 2 ) + D(x,w;a) 

lbn 2 



1 



16?r 2 
1 



1 



4tt 2 V 16tt 2 



V(G(w, 0; a) In G(w, 0; a)/M 2 ) + D(w, 0; a) 
V(G(x, 0; a) In G(a;, 0; a)/M 2 ) + D(x, 0; a) 



(1.13) 



On the other hand, the diagram in Fig. 3 displays a new kind of singularity, a generic 3-point one. The bare form 
of the amplitude reads 



G(x,0;a)G(y,0;a) (G(x,y;a)) 2 



(1.14) 



Let us recall Q that ordinary flat space renormalization goes in two steps, first renormalizing the inner 2-point 
singularity 



1 1 1 



111 \n(x-y) 2 M 2 



x 2 y 2 (x — y) 4 4 x 2 y 2 [x — y) 2 

Then, pulling derivatives out with exact manipulations, one finds 



(1.15) 



1 d I 1 d ln(x-y) 2 A/ 2 
4<9y M I x 2 y 2 dy^ (x - y) 2 



7r 2 ^)(y) 



lnx 2 M 2 



(1.16) 



This final piece needs the following supplementary, genuine 2-loop renormalization 



lnx 2 M 2 



2, ,2 



1 In x 2 M 2 + 2lnx 2 u 
□ o 



x 8 x' 

The coefficient in front of this last expression is responsible for yielding the correct 2-loop /3-function, 

3A 2 17 A 3 



/8(A) 



16tt 2 3 256tt 4 



(1.17) 



(1.18) 



once the first term of the anomalous dimension of 4> is introduced in the renormalization group equation. Coming 
back to our general discussion, it is easy to guess the correct generic renormalization of massless deformations of the 
flat theory. The first step is similar to the one in Eq. ( 1.15| ), 



G(x,0;a)G(y,0;a) 



167T 



■ V (G(x, y; a) In G(x, y; a)/M 2 ) + D(x, y; a 



(1.19) 



Through exact manipulations of covariant derivatives (recall that typically, T> = —g^V^d^ — a 2 , where a 2 is a 
constant), the expression will be transformed into a total derivative plus a generic 2-loop singularity which will be 
renormalized by 



G 2 {x) \nG{x)/M 2 



1 



32tt 2 



■V 



2 In 



G(x) 
M 2 



+ D'{x;M) 



D'(x; M) is again some left-over distribution. By taking the logarithmic M derivative of the above, we find 

M^D'(x;M) = -2D(x) , 



(1.20) 



(1.21) 



where D(x) is the distribution appearing in Eq. (1.8). This integrates to 
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D'(x; M) = D(x) \nG{x)/M 2 + f(x) , 



(1.22) 



where f(x) is some distribution independent of M. It turns out that this dependence is precisely what is needed for 
consistency of the RG equations for I^ 4 ) to order A 3 . This is another instance of the generality of our procedure. 
Both D(x) and f(x) give rise to benign singularities of the type 

1 1 1 (1.23) 



x 2 y 2 (x — y) 2 



which are integrable. Another fine point in Eq. ( 1.2C ) is that the relative coefficient —2 between the two powers of 
log is fixed by the identity itself and eventually produces the same 2-loop /3-function as in the flat case. The fact that 
we have renormalized logarithmic singularities has left no room for subleading divergences which might depend on a, 
leaving the result for the /3-function unaltered. It is, indeed, remarkable how universal these manipulations are. They 
substantiate the common lore about independence of some short-distance singularities on long-distance physics, while 
keeping closed explicit expressions for the whole amplitudes. 




Figure 4- 

Let us add a final and nontrivial example of this generic renormalization. The first non-local wave-function renor- 
malization can again be treated in full generality. The relevant diagram is shown in Fig. 4 and we shall call it the 
"setting sun" . Its generic renormalization is 

G 3 (x) = — —VCD + b)G( x) In G(x)/M 2 +D"(x). (1.24) 

5127T 4 

where 6 is a constant that depends on the details of the problem. It accompanies the term with a single T> op- 
erator, which is needed to cure subleading divergences which, for the first time, may introduce temperature and 
curvature dependences. As in all previous cases, D"{x) is a left-over distribution which carries information on the 
long-distance properties of the theory but does not affect short-distance renormalizations. A small effort applying the 
renormalization group equation yields the result for the leading anomalous dimension of the field <fi 

^>=M^r (L25) 

This generic renormalization again fulfills all three properties listed above: removal of singularities, correct renormal- 
ization group behavior and unitarity. 

Let us briefly discuss our ren orm alization procedure in the presence of masses. First of all, the equation for the free 
propagator changes from Eq. ( [L2] ) to 

(V + m 2 )G{x, 0; a; m) = "—^1 . (1.26) 



From this equation and the general form of the laplacian as previewed in Eq. (1.3), it becomes clear that the mass 
and the a parameter play similar roles. Their presence does not interfere with the leading singularities of the theory 
M. They first enter actively the g ame of renormalization to control the subleading renormalization of the 2-point 



function, which we called b in Eq. ( 1.24 ). Nonetheless, a blind application of massless recipes leads to unnecessarily 
complicated renormalized expressions due to the following observation. At long distances, the correct behavior of, 
say, the bubble diagram will be an exponential fall off as exp(— 2mr) since two particles circulate through the loop 
diagram. This is not the behavior of V(G log G/M 2 ), which decays as exp(— mr). Were we to use Eq. ([T^), the 
left-over distribution D would have to correct for this long-distance behavior. Although renormalization would be 
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carried out correctly, the answer achieved would not be minimally representing the long-distance fall-off. The solution 
to the problem is to substitute logG/Af 2 by a function which shares the same short-distance behaviour but provides 
a better representation of long-distance physics. A first step in this direction is to notice through examples that the 
more appropriate differential operator to be used to renormalize the bubble diagram turns out to be T> + Am 2 . A 
production threshold naturally finds its place next to the laplacian (and eventually sneaks in in momentum space 
as p 2 + Am 2 ), although it multiplies a distribution and thus there was no absolute need for it. The change in the 
distribution following the operator is of a more refined nature, and we pospone its discussion to the next section after 
an explicit example. Let us mention here that unitarity by itself demands further changes besides the one in the 
differential operator. 

We have just sketched a program to get explicit Green functions at low orders of perturbation theory for A</> 4 on 
modified spaces which leave the leading behavior of the propagator unchanged. The universality of the first two 
coefficients of the /3-function and the first one in the anomalous dimension of the field 4> also follows from the above. 
We shall now consider the explicit realization of these ideas. In Sec. 2, we treat constant curvature backgrounds, 
where we can in fact develop DR in an exact way and maintain manifest coordinate invariance, and in Sec. 3 we 
consider finite temperatures. 



II. RENORMALIZATION IN MAXIMALLY SYMMETRIC SPACES 



The renormalization effort in curved spaces, to date, has centered mainly on the techniques of point-splitting [Q, 
background field and heat kernel methods (||| • As opposed to flat space, in curved manifolds coordinate- rather than 
momentum-space is the natural setting for calculations, due to the lack of a generically-defined Fourier transform HQ. 
Because of that, it is reasonable to expect that differential renormalization will also be a natural method, beyond the 
abovementioned ones, since it is also based on coordinate space. The hope is that DR will present itself as an efficient 
renormalization procedure in curved space, and as we will see, this hope is brought to fruition. We also find that 
DR manifestly preserves the diffeomorphism invariance of amplitudes, given that all the typical manipulations with 
differential operators and their integration by parts are performed in this case in a covariant way. In this section, wc 
will concentrate on the simplest cases, namely, maximally symmetric spaces, and we will give here the renormalized 
expressions for a number of loop diagrams. In these spaces, the £i?<^ 2 coupling can be reabsorbed into the mass term 
and, thus, there is no need for a separate renormalization of £. 

a) The hyperboloid H n 

H n is the maximally symmetric Euclidean space of constant negative curvature in n dimensions. The simplest 
construction of maximally symmetric spaces is made through the embedding into an (n + I)-dimensional flat space, 
where the isometries of these spaces are very clearly put in evidence. Specifically for H n , this embedding is given by: 

^B2/V = -(y°) 2 + y / y i = -i, (2.f) 

er 

where y A ,A = 0, ...,n are coordinates in lR n+1 with metric t\ab — ( — V ■ ■■+), = 1, ...,n, and R = —n(n — l)a 2 is 
the (constant) Ricci scalar of H n . With this embedding, the metric in these spaces is given by 

ds 2 = (6» u -^f^)dy»dy\ (2.2) 
I + a z p z 

with y^ = 5^ v y v , and p 2 = y^y^. The action for a conformally coupled free scalar field is 

S = \j ' d n y \fg{y) (V(2/) d^(y) d^{y) + m 2 <j)(y) 2 + ^— ^ R 0(y) 2 ) . (2.3) 

Denoting by □ = g^ v V ^d u = — a 2 A the scalar Laplace-Beltrami operator, the equation satisfied by the scalar 
propagator in this theory is: 



2 n(n — 2)a 



-□ + (m z - ) 



G m (y) = -M- (2-4) 

V9(y) 



G 



In maximally symmetric spaces it is also possible to arrive at these propagators via a purely group-theoretical 
description based on highest weight representations of the isometry groups in these spaces ||7| . For the hyperboloid 
H n this group is SO(l, n), and if Lab, A,B = 0,1,..., n, are the corresponding generators, then highest weight 
representations are labeled by the eigenvalue 77 of L^ n (and possibly other eigenvalues related to the spin of the 
representation). Unitarity bounds 77 to be or > (n — 3)/2 in a continuous range. Invariance under boosts fixes 
uniquely the propagator, without reference to the free action for the corresponding field. However, we do not need the 
details of such a construction; rather, for our purposes here we just use the fact that 77 labels the scalar propagator 
just as the mass does, and we shall use this notation here for reasons which will become evident when we consider 
the renormalization of different diagrams. For now, the only difference that entails is that the equation satisfied by 
the scalar propagator, now labeled G n , reads 



[A + 77(77 - n + l)]G v {y) 



(2.5) 



where a 2 r]{rj — n + 1) = m 2 — n(n — 2)a 2 /4 is the quadratic Casimir eigenvalue for the representation 77. In order to 
find the propagator explicitly, we first write down the laplacian: 



a 2 A = (l + a 2 P 2 )-^ + 



(n-1) 



na p 



dp 



(2.6) 



It is possible to put this in a more workable form by defining the reduced variable z = y/l + a 2 p 2 , with range 



1 < z < 00. With this definition, Eq. (2.5) becomes 



i* 2 



s d 2 d 

1 )-Q^- nZ —+ f l(V 



dz 



1) 



G n (z) 



(2.7) 



where we write 8^ n '(y) in the original variables y for compactness, and where we have used \J g(y) \ y =o— 1. It is then 
straightforward to find that 

G n (z) = (z 2 -l)-^ 2 h v (z), (2.8) 



where h v satisfies an associated Legendre equation, with solutions given by Py(z) and Q%(z), the associated Legendre 
functions, with v = r\ — n/2 and p = (n — 2)/2 (for this and all further references to these functions, cf. ||). By 
examining the short distance limit z 1 of these functions it is easy to discard Pjf(z) as a solution since that leads 
to a propagator which is regular at the coincident point, and we finally find that the appropriate propagator is given 
by|: 



G v (z) 



2tt«/2 \2e™y/i^l 



(n-2)/2 



q { :~T(z)- 



(2.9) 



From now on we consider n = 4 for simplicity, although it is clear that analogous procedures will be possible in any 
number of dimensions. The normalized propagator becomes 



G v (z) 



1 



4^ 2 yfz 



(2.10) 



and it satisfies the equation 



77(77 -3)]G v (z) 



(2-11) 



In particular, the massless propagator obtains for 77 = 2 (and for 77 = 1, but this cannot be reached continuously from 
a massive theory, and we discard it here), and the equation it satisfies is: 



3 To be sure, this argument cannot discard a combination of both Q£ 
does. 



and P£ , but the highest-weight construction indeed 
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(2.12) 



We now proceed to find the DR identity appropriate to the diagram of Fig. 1. As in the case of massive fields in 
flat space ||, we first study these identities away from contact, and then include the appropriate contact terms which 
lead, in this case, to a well-defined flat limit. The simplest case, the massless bubble diagram of Fig. 1, exemplifies 
well the procedure: we attempt to write that amplitude as (covariant) differential operators acting on less singular 
quantities, and find the identity 



16ir 4 G 2 =2 (z) = 



a 2 zQl(z)Ql(z) 



z + 1 



(2.13) 



The differential operator on the r.h.s. is manifestly covariant, which implies we have preserved diffeomorphism 
invariance by this procedure. Also, for generic 77, Q] ) _ 2 {z) is expressed in terms of hyper geom etric functions which 
we cannot write in terms of simple functions. Thus, although we could have derived Eq. (2.13) more easily with the 
known explicit forms for Ql(z) and Qq(z), in general we cannot do that, and have to rely solely on recursion relations 
among these functions in order to find the appropriate DR identities. 

The contact terms for this DR identity will bring in a renormalization mass scale, and will determine the correct 
renormalization prescription for this diagram at all z > 1. They are gotten by examining the flat limit above 
{z — cosh ar, ap = sinh ar, where r is the geodesic distance to the origin in H n ): 



l^G 2 n=2 {z) 
a 2 Ql{z) 

zQ° (z) 
□ 



P 

z . z + 1 

— In - 

2 z- 1 



In a r 



2^2 



d z p + -d p + a\p z d 2 + A P d p ) 



□ flat 



leading to: 



1 

--□flat 



In a 2 r 2 



0. 



(2.14) 



(2.15) 



The r.h.s. of this is singular as a 2 — > 0, and we cure that by introducing a finite local quantity to the DR identity Eq. 



l^Gl =2 {z) = -(A -2) 



a 2 zQl{z)Ql{z) 



Vz^l 



M 2 

n 2 \n \ S^(y). 



(2.16) 



This then accomplishes the renormalization of the massless bubble diagram, with the usual DR prescription that the 
identity is now to be understood for all z > 1, and A as acting through integration by parts. The full, renormalized 
1-loop 4-point amplitude in Xcj) 4 theory is obtained by adding the tree- level amplitude to A 2 /2 times the above one, for 
each of s, t, u channels. Whe n thi s is required to satisfy the appropriate renormalization group equation, the correct 
value of the /3-function, Eq. (1.11), is obtained. Again, the simplicity of H n allowed us to circumvent an independent 
^-renormalization. 

At this point, this renormalized amplitude should be compared with the general procedure outlined in the previous 
section. There we exploited the fact that, generically, the short-distance behavior is the same as in flat space, and 
thus ultraviolet renormalization is achieved by "imi tating" as much as possible, given only G and T> generically, the 
renormalizati on o f the flat bubble diagram, Eq. (1/7). Thus, 1/x 2 is substituted by G, lnx 2 M 2 by — \nG/M 2 , and the 
result is Eq. fll.8| ) . That does accomplish the renormalization of G 2 generically, but at the cost of introducing an extra 
distribution D which corrects for the long-distance discrepancies this "imitation" process entails. Indeed, we have 
inde pende ntly checked that the same generic procedure works here, and does lead to the correct /3-function. However, 
Eq. (2.16) shows clearly that in this specific case we can do better than that, i.e., we have a specific, explicit expression 
valid at all distances without complicated distributions left over. This will not happen in the finite temperature case, 
as we shall see in Sec. 3, but does also happen for massive fields in flat space ||. The reason for this is that, loosely 
speaking, both in H4 and for massive fields, just like Q\ and K\ (the modified Bcssel function) are the respective 
equivalents of 1/x in flat space, so are Q° and K the respective equivalents of the flat space quantity In x 2 AI 2 , and 
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thus in these settings renormalization logarithms enter naturally while still allowing us to exploit recursion relations 
among these special functions. On the other hand, in these cases the renormalization scale will always enter separ ately , 
determined through a limi ting procedure as above. As a final comment, we note also that while M 2 enters Eq. ( 2.16| ) 
and the generic Eq. (1.6) in different ways, it does so in such a way that the logarithmic mass derivative of the 
renormalized amplitude gives a purely local term, which is uniquely responsible for the value of the 1-loop /^-function. 

We can now generalize the above to all £-loop 2-point functions of the form (G v —2(z)) +1 . It is not difficult to verify 
the following identities for the corresponding 2-, 3-, and 4- loop diagrams, for z ^ 1: 



——:Qo(z)) = -^(A + 4)(A-2)[a 2 z(Q 1 ,(z)) 2 Q°(z)] 

=QoW) = ^(A + 18)(A + 4)(A-2)[a 2 z V / ^T(QjW) 3 QS!W 



Vz 2 



?o(*) = 



18432 



(A + 40) (A + 18)(A + 4)(A - 2) [a 2 z(z 2 - l)(Qj(z)) 4 Q°(z)] 



(2.17) 



The first identity above, when supplemented by the appropriate contact term, represents the renormalized setting sun 
amplitude in X<p 4 theory. Also in A</> 4 theory, the second identity is a vacuum energy diagram, while the third one, 
being made of 5-point vertices, would not occur. For the setting sun diagram, the same discussion as above f or th e 
bubble diagram and its generic renormalization applies here , wit h the generic renormalization given in Eq. ( 1.24 ), 
and the ensuing anomalous dimension for <f> is given by Eq. (1.25). 

A pattern is now seen to emerge here which is valid for all t: for the ^-loop diagram, the differential operators which 
appear in the renormalization are given by: 



H[A + k v (k v -3)} 



(2.18) 



fc=i 



with r\ = 2, and the renormalizing quantity on the r.h.s. at each loop is written as 



a 2 z{z 2 -l)^r(Ql{z)YQl{z), 



(2.19) 



which in fact is always the same as the 1-loop one, from Qq(z) = — 1/vz 2 — T, but is written thus in order to preserve 
the same power of Legendre functions appearing on the l.h.s.. A comment is in order regarding the differential 
operators renormalizing these diagrams. In the fiat massive case ^], it was noted that these operators indicate the 
particle production thresholds for the corresponding diagrams, thus (□ — 4m 2 ) for the bubble, (□ — 9m 2 ) for the 
setting sun diagram, etc.. Here, a new but related phenomenon is seen to occur: if we consider the propagator s here 
to carry the representation rj = 2 (of the isometry group 50(1,4) of H4), then the differential operators in Eq. ( 2.18| ) 
will count the number of representations propagating along the diagram minus one. In curved spaces we cannot attach 
this the same meaning as in flat space, i.e., that this is an indication of the analiticity structure of the respective 
p-space amplitudes, but we nonetheless find this curious phenomenon worthy of mention. For conciseness, we do 
not present here the contact terms appropriate for the renormalization of these higher loop diagrams, but they are 
straightforward to find in the same fashion as for the bubble diagram. 

We can now move on to the renormalization of massive diagrams. They are somewhat more difficult to work out, 
and we present here the results only for the bubble and setting sun diagrams. In what follows, we need to use two 
recursion relations for associated Legendre functions H : 



dz 



1 



QZ+ 2 (z) =-2(// + l) 



41 2 



1 



(2.20) 



For the bubble diagram, we find after some manipulations: 



Vz 2 



t^)) 2 4( a+ 4- 2 



a 2 zQ\_ 2 {z)Q%_ 2 {z) 



Vz 2 ^! 



a 2 m 2 {Ql_ 2 {z)) 2 



(2.21) 



where m 2 = (77— 1)(?7 — 2)a 2 (compare Eq.(2.4) with (2.5) for n = 4 to find this). We note that, again, in the differential 
operator appearing here, although we cannot speak of an analiticity structure for the corresponding p-space amplitude, 
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we also see a "production threshold" for two particles of mass m appearing besides the curvature effect in the term 
—2 analyzed previously. Finally, to obtain the renormalized amplitude corresponding to this diagram, we need the 
appropriate contact terms. We find these in a two-step process: first, we take the limit a 2 — > while holding m 2 
fixed. This is done with the help of the following limit [10]: 



lim ^e-^Q^cosh £) = K^) , £ > 0, 

v— >oo 1/ 



where if„ are modified Bessel functions. By using this in the DR identity above, with z — cosh ar and v 2 
we find precisely the DR identity for massive scalars on flat space [||: 



/ mKi(mr) 
\ r 



1 . 9 . mKn(mr)Ki(mr) 
= - (D - 4m 2 ) ^ 



(2.22) 

m 2 /a 2 , 

(2.23) 



So far, no contact terms need be added. Now, we determine contact terms by the massless limit of the above equation. 
We simply borrow the result from Q, and finally write for the renormalized amplitude: 



VI 2 



!r,-2\ 



(*) = ^r A + 4— -2 



a 2 zgi_ 2 (z)Q°_ 2 (z) 



+ a 2 m 2 (Q° 2 (z)) 2 + tt 2 ln(16M 2 / 7 2 m 2 ) <5( 4 )(y) 



where 7 is the Euler-Mascheroni constant. Naturally, we could also have taken the direct a 2 — ► 0, m 2 



(2.24) 
limit to 



determine these contact terms, and the result would be identical to what we found for the massless case, Eq. ( 2.16 ). 
These two procedures differ by a choice of scheme or, more specifically, by the finite contact term 
7r 2 ln[(?7- l)(rj-2)j 2 /16]6^(y). 

We now present the DR identity that renormalizes the massive setting sun diagram. It represents somewhat more 
effort than previous cases, and we find: 



-16 



x/I 2- 



m 

A + ^-2 



{a 2 zQ%_ 2 {z){Q\_^z)f + m 2 z(Q°^ 2 (z)) 3 ) - 32a 4 mVz 2 - 1 (Q^z))' 



-8mV(a 2 + 8m 2 )y / z 2 ^l(Q^ 2 (z)) Q\_ 2 (z) - 28m 2 a 2 (a 2 zQ° v _ 2 (z)(Ql^ 2 (z)) 2 + m 2 z(Q°_ 2 (z)) 3 ) (2.25) 

By using the limit Eq. ( 2.22| ), the above identity again goes very smoothly in the a 2 — * 0, m 2 fixed limit into 
the corresponding amplitude in g|. Borrowing the massless limit of that expression, like for the setting sun above, 
determines the contact terms appropriate for this DR identity. The result above stands as a success of DR, insofar as 
it produces the fully renormalized 2-loop 2-point function in a compact and explicit way, something which had not 
been achieved previously. 

b) The sphere S n 



S n is the maximally symmetric space of constant positive curvature in n dimensions. Because it is a more familiar 
space, and because most formulas here would be essentially redundant with formulas from the previous section on the 
hyperboloid, we simply give here a more condensed summary of our results. Essentially all the geometry formulas 
from last section can be adjusted to spheres by taking a 2 —> —a 2 . The reduced variable to be used here, rather than z, 
is x = ±yl — a 2 p 2 , with range —1 < x < 1, where the positive(negative) sign holds for the upper(lower) hemisphere. 
The massive propagator is 



2tt"/2 \2e i7T VT 



(n-2)/2 



Q 



(n-2)/2 



(2.26) 



where now v{y + 1) = — m 2 /a 2 , and it satisfies the equation 



-(1 



dx 2 



d_ 

dx 



m 2 n(n — 2) 



G m (x) 



(2.27) 



10 



For the massless propagator in n — 4 in particular, we have: 

,2 j_ 



G m = (x) 



a 



4^ vr 



satisfying 



[A + 2] G m=0 (x) 



(2.28) 



(2.29) 



The renormalization of £-loop 2-point functions, analogously to Eqs. (2.13) and (2.17), works as follows: 



Vl — x- 

a 2 
Vl — x- 



y(A + 2) 



a 2 xQl{x)Ql{x) 
VT^x 2 



:Qh(x)) = 



--(A-4)(A + 2) [a?x{Qltx)fQl{xj\ 



Vl — x- 



Vl — x- 



- ; Ql(x)) = — (A -18)(A-4)(A + 2) [a^xVT^iQKx^Q ^) 



?Qo(*) = 



18432 



(A - 40) (A - 18)(A - 4)(A + 2) [a 2 a;(l - x 2 ){Ql(x)fQ°{x)} 



(2.30) 



The same pattern has emerged here as previously, with the only difference that all signs are reversed. Again, we 
do not work out the contact terms for these DR identities, as the procedure is identical to what was done previously. 
Also, one can verify that the first two DR identities above, which lead to the renormalization of the bubble and setting 
sun diagrams, respectively, again lead to the correct values for the /3-function and anomalous dimension in A</> 4 theory 
at lowest nontrivial order. 

With the hindsight of our work on the hypcrboloid, it becomes straightforward to check the renormalization identities 
for the massive bubble and setting sun amplitudes on the sphere. They are given, respectively, by: 



vr 



Ql(x)) =- A + 4 



a 2 xQl{x)Ql{x) 



vr 



a z m z {Ql(x)Y 



(2.31) 



and 



-16 



Vl — x 



- T Ql(x) 



r 2 

m 




A^ 


2 

m 




A + 9^- -4 






V2 


or 






'H 2 ' 





(a 2 xQl(x)(Ql{x)) 2 + m 2 x{Ql{x)f) + 32a i m 2 y/l - x* (Ql(x)f 
-8a 2 m 2 (a 2 - 8m 2 )Vl-x 2 (Q° u {x)) 2 Q%x) + 28m 2 a 2 [a 2 xQl (x)(Ql(x)) 2 + m 2 x(Q° l/ (x)) 3 ^ , 

and appropriate contact terms are again gotten by a straightforward massless and flat limit above. 



(2.32) 



III. RENORMALIZATION AT FINITE TEMPERATURE 



Although quantum field theory is usually defined at zero temperature, it is also possible to consider it at arbitrary 
nonzero temperatures |p], ll| , |l2| , and finite temperature effects are relevant, for instance, in studying phase transitions 
in the early universe |1S ] or the thermal structure of QCD |bj]. This implementation of temperature in quantum 
field theory is essentially made in one of two ways: either one proceeds in real Minskowskian time, or goes into 
Wick-rotated, compactified Euclidean time. The former suffers from some ambiguities |bj], while the latter is a more 
appropriate setting for our approach, and is the one we shall follow here. 

Again, to illustrate our general procedure, we work here with massless \(f> theory. It is known that finite tempera- 
tures induce masses through renormalization, and thus a consistent, renormalizable scalar model in this context must 
be massive. This phenomenon will be seen very clearly in our treatment from the non-closure of the renormalization 
group equations for the 2-point function. We choose, nonetheless, to start with a massless theory here in order to 
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illustrate our method in a simple way. Sure enough, the renormalization group equation will not be verified due to 
the lack of a bare mass to cancel a temperature-induced mass counterterm appearing at two loops. We then switch 
on a bare mass and treat it as a correction to the massless case in order to get its appropriate renormalization. The 
spirit of our approach remains, to show that in any theory one is able to isolate singular terms and correct them using 
differential identities. 

In imaginary time formalism, time is rotated into Euclidean time r and compactified to a cylinder of perimeter 

f>=±, (3-D 

where k stands for the Boltzmann constant. The momentum space propagator is 

1 2mi 

^2 ' Wn = 1 ~n~ > ( 3 - 2 ) 

w n -p P 

where the Euclidean zeroth component of the momentum vector takes only discrete values. In coordinate space this 
propagator can be written using the method of images as 

1 oo l 

G (^) = 4^ E x 2 + {T + n p )2 . ( 3 - 3 ) 

n— — oo v ! 

where x = Vx* 2 . This expression can be resummed, or, alternatively, the Fourier transform of the momentum 
propagator can be computed to yield 



G(x,r;(3)-- 



1 sinh 2£-x 



4nf3x cosh ^-x - cos 

v sinh vx 

8n 2 x cosh vx — cos vt 



(3.4) 



where we have rewritten the last line using v = =£■ = l^kT . Note that this compact form of the propagator is 
manifestly periodic in Euclidean time and decreases at long distances as 1/x, which is the characteristic decay in 
three dimensions. This is a hint at the effective reduction from four dimensions to three at high temperatures. 

The simplest diagram is a pure tadpole. Its departure from [3 = oo can easily be computed in coordinate space 



using the first form of the propagator, Eq. (3.3). The computation is reduced to 



lim (G(i,^-G(j,^ = M))=X;n 7 4 = i • (3.5) 

(x,t)-»(o,o) ^ o 4-7r z (np) z Ylp i 

This result agrees with the standard literature |0,|l3|. It can also be obtained through a simple expansion in Eq. 
( |3.4[ ). Nontrivial renormalization first appears when trying to compute the vertex function at 1-loop order. We need 
to renormalize G 2 , the square of the propagator, which presents a logarithmic singularity. We proceed as in Sec. 1, 
by noting that for x =/= 

G * = ttt^O (G(x InG (x ) M 2 ) + ——^G(x) - __ '— , 3.6 

lD7r z v lD7r^ 1287r 4 a; J smnw (cosh vx — cos vt) 

where G{x) is short-hand for the propagator and the laplacian operator is 

□ = d T d T + \d x X 2 3 x . (3.7) 
X 

The renormalization of this diagram has been carried out according to the generic analysis presented in Sec. 1. The 
singularity has been corrected in the laplacian term, and the last two terms of the r.h.s. above correspond to the 
expected left-over distribution, D, which encodes extra long-distance information of the theory. One can easily check 
that this distribution D has a leading singularity of the type [x 2 + t 2 ) -1 , which is Fourier-transformable. 

Now it is easy to compute the first coefficient of the p function, using the renormalization group equation for the 4- 
point amputated vertex function. At this order, there is no contribution coming from the anomalous dimension of the 



field <fi , thus Eq. (1.11) is readily obtained and we explicitly see that, to this order, the short-distance renormalization 



of the theory has not been affected by the compactification of Euclidean time. This corroborates the fact that we are 
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dealing with a scheme-independent term of the /3-function. To compute this term of the /3-function, only the coefficient 
of the logarithm matters, so we do not need the actual form of D(x), provided it is a well-behaved distribution. Yet, 
the point is that we have explicitly obtained the finite parts of the diagram without much effort. 

The 2-loop correction to the vertex function consists of two diagrams, Figs. 2,3. Fig. 2, as explained in Sec. 1, 
is just a convolution of the 1-loop diagram and only promotes logarithms to their squares, as dictated by the first 
coefficient of the /3-function, giving no contribution to the second coefficient of the /3-function. The o ther one is 
computed exactly as in Sec. 1, that is, in a first step we renormalize the inner div ergence th rough Eq. (1.19), with 
V = — □, and D given in Eq. (|3.6|) above. Identical manipulations lead to Eqs. (1.20,1.22), where now f(x) = 0. 



Again, the relative coefficient —2 between In and In above is dictated by short-distance behavior, and thus it is fixed 
and universal. Thus, collecting all the pieces, we finally obtain the renormalized expression 



yK-^^ + ii^ 



1 1 



(167T 2 ) 



rS (4) (x-y)D 



G(y) (ln^-21n^ 
! M 2 AP J 



(3.8) 



Where, again, we note that we have easily obtained explicitly the finite terms. 

We now turn to the self-energy correction depicted in Fig. 4, and again use the general recipe given in Sec. 1. 
Without undue effort, one finds 



G 3 (x) 



1 



512vT 4 



2v 2 )uG{x) In 



G{x) 
M 2 



R(x) , 



(3.9) 



where 



R(x) 



v I I (vx cosh vx — sinh vx) 
40967T 6 I x 3 sinh vx (cosh vx — cos vt) 



256tt 4 



G(x) 



1 (vx cosh vx — sinh vx) 2 2 sinh vx + vx cosh vx — vx cos vt 



2048^ 6 x 4 ( C osh vx — COS VT) 



sinh; 



(3.10) 



This has the same short-distance behavior as the propagator, and is Fourier-transformable. 

As explained in Sec. 1, the appearance of a term linear in the laplacian (with coefficient 2v 2 ) has to do with 
subleading singularities. In this case, it reflects the fact that the thermal bath screens the interaction. This will 
generate the inconsistency we are expecting: the renormalized 2-loop 2-point function will be 



r< 2 ) = -nS^(x) 



x 2 ( 1 



2v 2 )uG(x) In 



G(x) 



6 V512tt 4 V " ' " /""V"/™ M 2 

and to check whether it satisfies the appropriate RG equation we compute 

v2 / 1 X 2 



R(x) 



dM 6 V 16tt 2 



(p5 (4) (x) + 2v 2 S^(a 



(3.11) 



(3.12) 



The first term leads to the expected value of the anomalous dimension, Eq. (1.25), while the second term represents 
a mass renormalization. That means that, as mentioned previously, renormalizability requires a massive theory ab 
initio. 

It is clear that masses modify the long-distance behavior of the propagator while preserving the 1/x 2 singularity at 
the origin. Since renormalization consists in curing short-distance singularities, the possibility is open to treat masses 
as perturbations. More precisely, the exact massive propagator at finite temperature can be found to be 



G(x, t; /3; m) 



8n 2 x 



E 



^—invT — x\J v 2 n 2 +m 2 



(3.13) 



Expanding around the massless case, we find 



G(x, t; P; m) 



sinh vx 



8tt 2 x V cosh vx — cos vt 



2v 



In [v 2 



2 )] 



(3.14) 
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We are now ready to use this modified expression of the propagator to go through the renormalization procedure in 
all previous graphs. It is readily noticed that the 1- and 2-loop coefficient of the /3-function and the 2-loop coefficient 
of the anomalous dimension of the scalar field are not modified by the presence of masses. The first and relevant 
appearance of a mass countcrtcrm takes place in the setting sun diagram, as explained above. The expression for that 
diagram corresponds to Eq. ( |3. 11 ) with the addition of the first massive correction 



16tt 2 



, , G(x) m 2 
u mU G{x)hx^ r + — U 



G(x) In 2 [v 2 (x 2 + r 2 )] -21n 



G(x) 
M 2 



(3.15) 



The complete 2-point function then verifies the renormalization group equation 



M m + ^x + ^ m2 ^- 2 ~ 2 ^ r{2) ^-° 



(3.16) 



with 



7m 



16tt 5 



3 \mJ m b I 



(3.17) 



which is a scheme-dependent coefficient. 



IV. CONCLUSIONS 



The aim of the present paper has been to exploit the idea that the renormalization process amounts to replacing 
short-distance singularities by distributions and, more importantly, that this project can be carried out for any 
background a field theory might be defined on. Explicit examples on constant curvature and thermal backgrounds 
have been worked out. 

We can draw several lessons from this investigation. In the finite temperature case, the absence of a compact 
form for the massive propagator in position space has forced us to treat masses in a perturbative way. This is a 
drawback, though renormalization group constants are obtained probably more easily than with other methods. On 
the other hand, constant curvature backgrounds are surprisingly suited for differential renormalization. Explicit closed 
expressions for 2-loop amplitudes are obtained with ease. 

Further applications of this renormalization procedure are manifold, as for instance the extension to QED, along 
the lines of |lq] , or computations of the effective potential on curved backgrounds. Its very essence, based on the 
coordinate space propagator, makes it an extremely simple method to renormalize a field theory. 
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